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Abstract 

Using the 1.9 pb~ l of data collected with the CMD-2 detector at 
VEPP-2M the decay mode 4> ~^ Wli V ~^ 7r + 7r~7r° has been studied. 
The obtained branching ratio is B(</> — > 777) = (1.18 ± 0.03 db 0.06)%. 



1 Introduction 



Radiative transitions between vector and pseudoscalar mesons are very 
interesting for tests of the quark model, SU(3) symmetry and Vector Domi- 
nance Model 0. Recently there were extensive discussions concerning mech- 
anisms of SU(3) breaking, the role of anomalies and a possible admixture 
of glue in mesons in radiative decays [0, ^|, |], ^ ||, [j], ||. Despite numerous 
experimental efforts, many open questions still exist requiring more precise 
measurements of the decay probabilities. 

The radiative magnetic dipole transition of into rj has been previously 
studied in many experiments M based on the neutral final states arising 
when f] decays into 77 or 37r°. In this work, using the general purpose CMD- 
2 detector at the high luminosity Novosibirsk collider VEPP-2M, the <fi — > 777 
decay rate into the charged-neutral final state was measured: e + e~ — > <p — > 
777 — > 7r + 7r~7r°7. The study is based on a data sample of about 1.9 pb -1 
corresponding to 3.5 million meson decays collected in 1996 with CMD-2. 



2 Experiment 



The CMD-2 detector is described in detail elsewhere p0[ |TT| , [12|]. It is an 
axial field detector consisting of a drift chamber with about 250 /i resolution 
transverse to the beam and proportional Z-chamber used for trigger, both 
inside a thin (0.4 Xq) superconducting solenoid with a field of 1 T. 

The barrel calorimeter placed outside of the solenoid consists of 892 Csl 
crystals of 6 x 6 x 15 cm 3 size and covers polar angles from 46° to 132°. 
The energy resolution for photons in the Csl calorimeter is about 9% in the 
energy range from 50 to 600 MeV. 

The BGO (B^GeaO^) end-cap calorimeter consists of two identical end- 
caps and covers forward-backward angles from 16° to 49° and from 131° 
to 164°. It is placed inside the main solenoid in the gap between the drift 
chamber and the flux return yoke and operates in a 1 T longitudinal magnetic 
field. Each end-cap is a compact matrix of rectangular BGO crystals of 
2.5x2.5x15 cm 3 size. The thickness of the calorimeter for normally incident 
particles is equal to 13.4 X . The calorimeter consists of 680 BGO crystals 
with a total weight of about 450 kg. The light from each crystal is read 
out by vacuum phototriodes which can operate in magnetic fields up to 2 T. 
The energy and angular resolution of the end-cap calorimeter is equal to 
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a E /E = 4.6%/^fE{GeV) and a^ e = 2 • 10~ 2 / 1 E(GeV) radians respectively. 

Data was taken in the 0-meson energy range in April - July, 1996. It was 
the first run after installing the BGO end-cap calorimeter. The information 
about the data sample included in the analysis is presented in Table |I[ The 
luminosity was determined using the number of detected e + e~ — > e + e~ events 
0. 



3 Selection criteria 

Events with two tracks and more than one photon were selected using 
the following criteria: 

• One vertex is found in the event 

• Two tracks with opposite charges are reconstructed from this vertex 
and there are no other tracks 

• The angles of both tracks with respect to the beam are limited by 
40° < 9 < 140° to match the optimal drift chamber coverage 

• The number of photons detected in the Csl and BGO calorimeters is 
more than one and less than six. The cluster in the calorimeter was 
accepted as a photon when it did not match any charged track and its 
energy was more than 30 MeV in the Csl calorimeter (more than 40 
MeV in the BGO calorimeter). 

• The distance from each track to the beam R m i n < 0.2 cm 

• The distance from the vertex to the interaction point along the beam 
direction \Z vert \ < 10 cm 

• The space angle between the tracks Aip < 143° to suppress events of 

- K S K L , K S - 7T + 7T- 

• The angle between the tracks in the v-tp plane A(p < 172° 

• The total energy of the charged particles (assuming that both are 
charged pions) 6^+^- < 520 MeV 

• The energy deposition of the photon accepted as a recoil one is more 
than 250 MeV 
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The cut for the total energy 6^+^- < 520 MeV is to reject events from 
the decay (ft — *■ 7r + 7r~7r°. From simple kinematical relations the total energy 
of charged pions from the (ft — > 7r + 7r~7r° is greater than 539 MeV while for 
</> ^ ??7 it is less than 522 MeV. 

4 Kinematics of <f> — ► 777 

In this study the final state contains two charged particles and three 
photons. 

Since (ft —>■ 7/7 is a two-body decay and 77 is a narrow state, the energy of 
the recoil photon is given by 

(2^) 2 ~ m* 

^ = AE b ' (1) 

where is the beam energy, is the r\ mass. The invariant mass of two 
other final photons is equal to m^o. At the (ft meson peak the energy of the 
recoil photon equals 363 MeV. 

In this analysis detection of all final particles is not necessary. Instead, 
one requires the detection of both charged pions and at least one photon, the 
recoil one. Using the measured momenta of charged pions, angles <p and 9 of 
the recoil photon and the assumption that the energy of the recoil photon is 
given by ([I]), one can reconstruct the invariant mass of all other photons in 
the system, M inv . In case of the process e + e~ — ► (ft —* r]j,r] — ► 7r + 7r~7r this 
parameter M inv = m^o . 

The recoil photon has the largest energy among three photons in the final 
state. So, the simplest way would be to take a photon with the largest energy 
deposition in the calorimeter. In this study, however, we prefer not to use 
the energy deposition at all, relying instead on the position of the cluster 
(angles of photons). The Monte Carlo simulation of the decay mode studied 
in this work shows that the recoil photon prefers to fly opposite to the total 
momentum of charged pions p n + n - in contrast to the photons from the 7r° 
decay (Fig. [1]). Therefore the photon whose direction is closest to — p 7T + 7T - 
was selected as the recoil one. 
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cosa cosa 

Figure 1: Distributions over cosine of the photon angle with respect to 
—p w + n - for the simulation of <fi — > 777 — > 7r + 7r~7r°7: a) for the recoil pho- 
ton; b) for one of the photons from the n° decay. 
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Figure 2: Mf nv distribution for selected experimental data. 
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Figure 3: Mf nv distribution for data (points with errors) together with the 
sum of simulated background processes (hatched histogram). The contribu- 
tion of — > 7T + 7r~7r° is shown by the solid hatched histogram. 



5 Data Analysis 

The distribution over the parameter Mf nv for selected events is presented 
in Fig. |||. Events from the decay <fi — > 777 produce a peak near the ir° 
mass squared. Other events surviving after all the cuts are coming from the 
following processes: 

• <fi — > 777, 77 — ► 7T+7r - 7. For such events the value of M inv should be 
obviously zero and a corresponding peak can be seen in Fig. ^|. 

• e + e~ — > urn — ► 7r + 7r _ 7r°7r° 

• (f) — ► 7r + 7r _ 7r . Most of the events from this decay are suppressed by 
the cut on the total energy of the charged particles e n + 7T - < 520 MeV, 
but some small fraction could survive. 

• 4> K S K L 

• beam background, the photon conversion at the vacuum tube 

Simulation of the first three background processes used the GEANT based 
CMD-2 simulation code WM. It is difficult to study two last sources of 
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background using simulation, therefore their behaviour was studied with the 
help of experimental data using the selection R m i n > 0.3 cm instead of 
Rmin < 0.2 cm. Figure |] presents the distribution over Mf nv for simulated 
events of background processes after applying selection criteria. The fraction 
of each background process in the total number of simulated events corre- 
sponds to its cross section, the fraction of background with large distances 
from the beam has been also normalized. The sum of background distribu- 
tions is shown by the hatched histogram in Fig. 0. The logarithmic scale was 
chosen for better demonstration of the background. The background simu- 
lation is in reasonable agreement with the data supporting our description 
of the background processes, sources is correct. The background from the 
— > 7r + 7r~7T decay (see the solid hatched histogram in Fig. ^ is rather small 
but unpleasant since its distribution for events surviving the selection has a 
maximum near that of the signal. 

The number of — > 777 events in the mode 77 — > 7r + 7r _ 7r was determined 
by the following procedure: 

• The distribution over Mf nv for each background process was fit by a 
smooth function to fix its shape. 

• The distribution over Mf nv for all the data (Fig. 0) was fit by the sum 
of background functions excluding the signal region. 

• After that the function describing the shape of the signal was deter- 
mined by the fit including the signal region and taking the sum of 
background functions with fixed parameters. 

• As a result of previous steps one had the functions describing the shape 
of signal and each background component. For each beam energy the 
number of — > 777 — ► 7r + 7r~7r°7 events was determined by fitting the 
distribution over Mf nv by the sum of these functions with the fixed 
shape and the number of events for signal and background processes 
as free parameters. For the — > 7r + 7r~7r° background the detection 
efficiency was determined from the simulation and the number of events 
at each energy was calculated according to the energy dependence of 
the cross section of the reaction e + e~ — > — > tt + 7t~7t° jHJ. The total 
background of the decay mode — > tt + 7t~7t° is less than 2%. 

The resulting number of 777 events and the obtained cross section are 
presented in Table [I]. 
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At each energy the cross section was calculated as 

N 
1 '777 

Q = — 

C ■ (1 + 5) • £ ■ Bjj^+ir-Tr ' 

where C is the integrated luminosity at this energy, 5 is a radiative correction 
for the process — > 777, e is the detection efficiency, B v ^ n + w - n o = 0.231 ± 
0.005 §. 

The radiative correction was calculated according to [16| and its behaviour 
in the 0-meson energy range is presented in Fig. f|. 




Figure 4: Radiative correction 5 for <p — > 777 vs energy 

The detection efficiency was determined using simulation of — > 777, r] — > 
7t + 7t~7t° and the value e = 0.247±0. 001 ±0.007 was obtained. The systematic 
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uncertainty of the efficiency was estimated by comparing results of the simu- 
lation of nuclear interactions by different methods (FLUKA and GHEISHA) 
[T4|. For simulation of the decay r) — > tt + 7t~7t° instead of standard GEANT 
tools which assume a constant matrix element we used the following matrix 
element parametrisation from [ 17|j 



\M(x, y)\ 2 oc (1 + ay + by 2 + cx + dx 2 + exy). 
Here x = VS(T + -T_)/Q,y= (3T /Q) - 1, 

where T+, T_ and To are the kinetic energies of the tt + , 7r~ and 7r° in the r\ 
rest frame, and Q = T + + T_ + Tq. The parameter values were taken from 



recent Crystal Barrel data [T8 



The cross section of the process e + e~ — > <p — * VI was parametrized by a 
sum of the 0-meson amplitude and an interfering term with the tails of the 
p and to resonances: 



a 



VI ~ s 3/2 



\^6 + A\ 2 , 



A, 



s-m 2 +iyfsT ^(s) ' 



where s = 4E 2 , m<f>, T^, cr2 are the mass, width and peak cross section 
(s = for the 0-meson, and F rry (s) = (y/s(l — m 2 /s)) 3 is a function 
describing the dynamics of the (j) — *■ 777 decay including the phase space. 

The energy dependence of the <fr meson width was written as 



r*(s) = • (b k+k - + B 



mlF K+K - (s) m 2 F KsKL 



\ A ' A s- F k +k- (mj) KsKh s ■ F KsKl (m 
m^F^-^m^) 

where B K + K -, Bk s k l , B^+^-^o and B^, are branching ratios of the me- 
son decay into corresponding modes; F K ^(s) = (s/4 — m^) 3 / 2 , and for the 
F^+ n -^o (s) calculation the model assuming the decay </>(c<j) — > p7r — > 7r + 7r~7r° 
was used. 

Fitting the energy dependence of the cross section of — > 777 (Fig. |5|) the 
meson mass and peak cross section were determined, while was fixed 
at the world average value || . The following results were obtained when the 
interfering term A = 0: 

erg = 49.8±1.2±2.4 nb; m = 1019.38±0.07±0.08 MeV; x 2 /U/- = 18.4/16. 
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Table 1: The c.m.energy, integrated luminosity, number of 777 events, radia- 
tive correction and cross section of — > 777 



No. 


Vs, MeV 


JCdt, nb' 1 


N 


8 


(Tyyy, nb 


1 


985.80 ±0.20 


93.4 ±0.5 


3.2 ±3.8 


-0.21 


0.8 ±0.9 


2 


1004.80 ±0.20 


98.8 ±0.5 


7.6 ±4.3 


-0.24 


1.8 ± 1.0 


3 


1011.44 ±0.20 


108.9 ±0.6 


11.8 ±4.7 


-0.25 


2.5 ± 1.0 


4 


1016.14 ±0.20 


123.4 ±0.6 


110.0 ± 11.4 


-0.28 


21.7 ±2.3 


5 


1017.08 ±0.06 


177.8 ±0.7 


192.1 ± 14.9 


-0.28 


26.3 ±2.0 


6 


1018.04 ±0.06 


226.2 ±0.8 


347.4 ± 20.0 


-0.28 


37.4 ±2.2 


7 


1018.88 ±0.04 


84.6 ±0.5 


170.3 ± 14.1 


-0.27 


48.4 ±4.0 


8 


1019.04 ±0.04 


63.8±0..4 


128.9 ± 12.3 


-0.27 


48.5 ±4.6 


9 


1019.20 ±0.04 


103.4 ±0.5 


239.9 ± 16.6 


-0.26 


55.0 ±3.8 


10 


1019.84 ±0.04 


87.0 ±0.5 


185.4 ± 14.5 


-0.24 


49.2 ±3.8 


11 


1020.08 ±0.04 


170.7 ±0.7 


340.5 ± 19.6 


-0.22 


44.9 ±2.6 


12 


1020.72 ±0.04 


40.8 ±0.3 


57.6 ±8.5 


-0.19 


30.6 ±4.5 


13 


1020.98 ±0.06 


143.8 ±0.6 


207.7 ± 15.7 


-0.17 


30.5 ±2.3 


14 


1021.80 ±0.06 


119.5 ±0.6 


102.1 ± 11.4 


-0.11 


16.8 ± 1.9 


15 


1022.72 ±0.06 


77.6 ±0.5 


56.4 ±8.3 


-0.04 


13.3 ±2.0 


16 


1028.26 ±0.10 


84.7 ±0.5 


9.2 ±4.0 


0.28 


1.5 ±0.6 


17 


1033.70 ±0.18 


55.7 ±0.5 


3.2 ±3.5 


-0.01 


1.0 ± 1.1 


18 


1038.96 ±0.12 


22.9 ±0.5 


2.5 ±3.5 


-0.09 


2.1 ±2.9 
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Figure 5: The cross section of e + e — > <\> — > 777 with the fit function. 

If A is calculated accurately taking into account the p and u> meson con- 
tributions, the value of <r° is changed by only 1%. 

The systematic uncertainty of the peak cross section is estimated to be 
4.8% and includes the following contributions: 

• luminosity determination (2%); 

• the uncertainty of the value of B(rj — > 7r + 7r~7r°) (2.2%); 

• determination of the detection efficiency (2.8%); 

• uncertainty of the background subtraction (2.5%) 

The systematic uncertainty of the mass was estimated to be 80 keV from 
the difference in the c.m.energy obtained from two different methods of the 
beam energy determination (one used the momenta of charged kaons and the 
other was based on the analysis of the magnetic field of the collider [|13|] ) . 
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Taking the value T ee /T tot = (2.99 ± 0.08) x 10~ 4 from §, one obtains 
the following value of the branching ratio corresponding to the peak cross 
section above: 

B^ rn = (1.18 ± 0.03 ± 0.06)%, 

where the first error is statistical and the second one is systematic. The 
latter includes contributions from the uncertainties of the peak cross section 
(4.8%) and that of T ee /T tot (2.7%). 

6 Conclusions 

Using a data sample of about 3.5 million meson decays the CMD-2 
collaboration studied the decay — > 777, 77 — > 7r + 7r~7r°. About 2200 events 
of this decay have been selected. The (f) meson mass was determined to be 
(1019.38 ± 0.07 ± 0.08)MeV and agrees with other measurements |§. The 
branching ratio of the decay was measured with high precision to be (1.18 ± 
0.03 ± 0.06)% and is consistent with the world average value (1.26 ± 0.06)% 
as well as with another accurate determination of this quantity recently 
reported by the SND collaboration (1.246 ± 0.025 ± 0.057)% @. It is the 
first measurement of the <fi — > 777 decay rate using the charged mode of the 
77 meson decay. Analysis of a five times bigger data sample is in progress. 
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